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We create molecules from fermionic atoms in a three-dimensional optical lattice using a Feshbach 
resonance. In the limit of low tunnelling, the individual wells can be regarded as independent three- 
dimensional harmonic oscillators. The measured binding energies for varying scattering length agree 
excellently with the theoretical prediction for two interacting atoms in a harmonic oscillator. We 
demonstrate that the formation of molecules can be used to measure the occupancy of the lattice 
and perform thermometry. 
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Quantum degenerate atomic gases trapped in the peri- 
odic potential of an optical lattice form a quantum many- 
body system of unprecedented purity. The short-range 
interaction due to atom-atom collisions makes optical 
lattices ideal to experimentally realize Hubbard models 
ll y] . Experimental studies of bosonic Mott insulators 
& DL 0- El an d of a fermionic band insulator \j] have 
provided a first taste of this new approach to quantum 
many-body physics. 

In the vicinity of a Feshbach resonance the collisional 
interaction strength between two atoms is tunable over 
a wide range. For two fermionic atoms on one lattice 
site strong interactions change the properties of the sys- 
tem qualitatively and physics beyond the standard Hub- 
bard model becomes accessible. Crossing the Feshbach 
resonance in one direction leads to an interaction in- 
duced coupling between Bloch bands which has been ob- 
served experimentally Q and described theoretically || . 
Crossing the resonance in the other direction converts 
fermionic atoms into bosonic molecules. These processes 
have no counterpart in standard condensed matter sys- 
tems and demand novel approaches to understand the 
mixed world of fermions and bosons in optical lattices 
E3 EH El • Descriptions based on multi-band Hub- 
bard models are extremely difficult to handle and there- 
fore the low-tunneling limit is often used as an approxi- 
mation. In this limit the lattice is considered as an array 
of microscopic harmonic traps each occupied with two 
interacting atoms in different spin states. 

The harmonic oscillator with two interacting atoms 
has been studied theoretically and the eigenenergies have 
been calculated in various approximations 0, Il3lll3 . Hif . 
Its physics is governed by several length scales. The 
shortest scale is the characteristic length of the van- 
der-Waals interaction potential between the atoms. The 
next larger length scale is given by the s-wave scatter- 
ing length characterizing low-energy atomic collisions. 
However, near the Feshbach resonance it may become 
much larger than the extension of the harmonic oscilla- 
tor ground state. A precise understanding of the interac- 
tions in this elementary model is a prerequisite in order to 
comprehend the many-body physics occurring in optical 



lattice systems with resonantly enhanced interactions. 

In this paper we study a spin-mixture of fermionic 
atoms in an optical lattice and their conversion into 
molecules by means of a Feshbach resonance. The bind- 
ing energy as a function of the s-wave scattering length 
between the particles is measured and compared with 
theoretical predictions. Moreover, we demonstrate that 
the molecule formation can serve as a measure of the tem- 
perature of the atoms in the lattice. For mapping out the 
phase diagram of many-body quantum states in the lat- 
tice the temperature is a key parameter. So far, temper- 
ature has not been measured in a lattice since standard 
methods - such as observing the rounding-off of the Fermi 
surface - turned out to be dominated by the inhomogene- 
ity of the trapping potential rather than by temperature 



However, we find that the occupancy of the lattice 
depends strongly on the temperature LLa, LUj and the 
conversion of pairs of atoms into molecules is a sensi- 
tive probe, similar to the case of harmonically trapped 
fermions |20j. 

In a previous experiment, deeply bound molecules of 
bosonic atoms have been created in an optical lattice by 
hoto-association and were detected by a loss of atoms 
*2l) . There the binding energy is only determined 
by the atomic properties and does not depend on the 
external potential, nor can the scattering properties be- 
tween the molecules be adjusted. In low-dimensional sys- 
tems molecules produced by a Feshbach resonance using 
fermionic atoms have been observed recently [22l . 

Our experimental procedure used to produce a degen- 
erate Fermi gas has been described in detail in previous 
work 0- In brief, fermionic 40 K atoms are sympatheti- 
cally cooled by thermal contact with bosonic 87 Rb atoms, 
the latter being subjected to forced microwave evapora- 
tion. After reaching quantum degeneracy for both species 
we remove all rubidium atoms from the trap. The potas- 
sium atoms are then transferred from the magnetic trap 
into a crossed-beam optical dipole trap which consists of 
two horizontally intersecting laser beams. In the optical 
trap we prepare a mixture of the \F — 9/2, mp = —9/2) 
and the \F = 9/2jn F = -7/2} state with (50 ± 4)% in 
each spin state |2 1| and perform additional evaporative 
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cooling at a magnetic field of B = 227 G. At the end of 
the evaporation we reach temperatures of T/Tp = 0.25 
with up to 2 x 10 5 particles. The temperatures are deter- 
mined from a fit to the density distribution of the non- 
interacting atomic cloud after ballistic expansion. 

We tune the magnetic field to B — (210.0±0.1) G, such 
that the s-wave scattering length between the two states 
vanishes. The atoms are then transferred into the optical 
lattice formed by three orthogonal standing waves with a 
laser wavelength of A = 826 nm Q. The resulting optical 
potential depth Vq is proportional to the laser intensity 
and is conveniently expressed in terms of the recoil energy 
E r = h 2 k 2 /(2m), with k — 2ir/\ and m being the atomic 
mass. 




Pulse detuning 8 [kHz] 
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FIG. 1: a) Illustration of the rf spectroscopy between two 
bound states within a single well of the optical lattice. The 
atoms in the initial states | — 7/2} and | — 9/2) are converted 
into a bound dimer by sweeping across a Feshbach resonance. 
Subsequently we drive an rf transition | — 7/2) — » | — 5/2) to 
dissociate the molecule, b) rf spectrum taken at B = 202.9 G, 
i.e. for a < 0. c) rf spectrum taken at B — 202.0 G, i.e. for 
a > 0. All data are taken for a lattice depth of 22_E r . The 
lines are Lorentzian fits to the data. 

In the optical lattice we create a band insulator for each 
of the two fermionic spin states 0. Subsequently, the 
molecules are formed by ramping the magnetic field from 
the zero crossing of the scattering length at B = 210 G in 
10 ms to its final value close to the Feshbach resonance 
located at Bq — 202.1 G [2^. From the parameters of 



our magnetic field ramp we estimate that the molecule 
formation is performed adiabatically. We measure the 
binding ener gy Ep o f the dimers by radio- frequency spec- 
troscopy |23l l26l l27| . A pulse with a frequency vrf and 
a duration of 40 fis dissociates the molecules and trans- 
fers atoms from the state | — 7/2) into the initially un- 
populated state | — 5/2) which does not exhibit a Fesh- 
bach resonance with the state | — 9/2) at this magnetic 
field. Therefore the fragments after dissociation are es- 
sentially noninteracting (see fig. QJi) . We vary the detun- 
ing 5 = vjif — vq from the resonance frequency i>o of the 
atomic | — 7/2) — > | — 5/2) transition. The power and du- 
ration of the pulse is chosen to constitute approximately 
a 7r-pulse on the free atom transition. The number of 
atoms in each spin state is determined using absorption 
imaging after ballistic expansion. For this we ramp down 
the lattice exponentially with a duration of 1 ms and time 
constant of 0.5 ms from the initial depth Vo to 5 E r to re- 
duce the kinetic energy of the gas and then quickly turn 
off the trapping potential within a few /zs. The magnetic 
offset field is switched off at the start of the expansion, so 
that no molecules can be formed in the short time that 
it passes the Feshbach resonance. We apply a magnetic 
field gradient during 3 ms of the total 7 ms of ballistic 
expansion to spatially separate the spin components. 

Figure^ and c show rf spectra of atoms and molecules 
trapped in a three-dimensional lattice with a potential 
depth of of Vq — 22 E r , which corresponds to a trapping 
frequency of uj = 2ir x 65 kHz in the potential wells of the 
lattice. The spectrum in figure ^ is taken at a magnetic 
field of B = 202.9 G, corresponding to a/ato = —1.3. We 
have calculated the ground state radius aho = 64 nm by 
minimizing the energy of a Gaussian trial wave function 
inside a single well of our lattice potential. As compared 
to a Taylor expansion of the sinusoidal potential around 
the minimum this results in slightly more accurate values. 
The s-wave scattering length is denoted by a. For nega- 
tive scattering length the molecules are only bound when 
they are strongly confined whereas no bound state would 
exist in the homogeneous case. The spectrum exhibits 
two resonances: the one at S = corresponds to the 
atomic transition from the | — 7/2) into the | — 5/2) state. 
This transition takes place at all lattice sites which ini- 
tially were only singly occupied and no molecule has been 
formed. The second resonance at <5 > corresponds to 
the molecular dissociation and is shifted from the atomic 
resonance by the binding energy. Simultaneous with the 
increase in the | — 5/2) atom number we observe a loss of 
atoms in the | — 7/2) state, whereas the | — 9/2) remains 
unaffected. This is expected since the rampdown of the 
lattice before detection dissociates all molecules and the 
| — 9/2) atom number should be fully recovered. Resid- 
ual fluctuations of the | — 9/2) atom number are probably 
due to uncertainties in the atom number determination 
since we do not observe a specific pattern in our spectra. 

With the magnetic field tuned to B — 202.0 G (see fig. 
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H 3 ), which corresponds to a/a^o = 16.8, the spectrum 
changes qualitatively. For this value of the scattering 
length stable molecules exist even in free space but the 
molecules formed in the lattice are not detected by our 
state-selective imaging procedure unless they are dissoci- 
ated by the rf-pulsc. Therefore only rf dissociated atom 
pairs show up in the time-of-flight images, resulting in 
an increasing number of atoms in the | — 5/2) and the 
| — 9 /2) state at the molecular resonance. 

In contrast to earlier work, where molecules were dis- 
sociated into a continuum and the fragments were es- 
sentially free particles |23l l26l |2jj , the fragments in our 
configuration occupy an energy eigenstate of the confin- 
ing potential. In such a bound-bound transition no extra 
kinetic energy is imparted onto the dissociated fragments 
since any excess excitation energy would have to match 
the band gap. We determine the binding energy from the 
separation of the atomic and the molecular peak. More- 
over, since there is at most one molecule present per lat- 
tice site, collisional shifts [2^, |2^| are absent and we can 
estimate the error in the binding energy from the fit error 
which is less than 5 kHz. 
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FIG. 2: The measured binding energy of molecules in a three- 
dimensional optical lattice. The data are taken for several po- 
tential depths of the optical lattice of 6 E r (triangles), 10i? r 
(stars), 15 E r (circles), and 22 E r (squares). The solid line 
corresponds to the theory of ref. [l^l with no free parameters, 
the dashed line uses an energy-dependent pseudopotential ac- 
cording to ref. [T^| . At the position of the Feshbach resonance 
(a — > ±00) the binding energy takes the value E = —Tiu. 

We have investigated the dependence of the binding 
energy of the molecules on the scattering length (Fig.[2J). 
The scattering length is derived from the magnetic field 
using the parametrization of the Feshbach resonance 
a(B) = dhnCL — b^Bo ) ' with a^g = 174 ao H3 and 
AB = 7.8 G |3J. We compare our data with the theory 
for two particles trapped in a harmonic oscillator poten- 
tial interacting via an energy-independent pseudopoten- 
tial [13 . The binding energy E of the molecules depends 



on the scattering length according to 

r- T(-E/2hu) a ho 
T(-E/2hu -1/2) a ' [ ' 

where T(x) denotes the Gamma function. We find the 
normalized binding energy Ejhio to be independent of 
the strength of the lattice and all data points agree well 
with the theoretical prediction of equation (1) without 
adjustable parameters. A pseudopotential approxima- 
tion is valid, as long as a; lo is large compared to the 
characteristic length scale of the van-der-Waals poten- 
tial between the two atoms (3§ = (mCe/h 2 ) 1 ^ [l4L fl5| . 
which for our experiments is a/j / (3$ > 10. However, for 
a 3> &ho an energy-dependent pseudopotential will model 
the system more accurately. We have calculated the ef- 
fective range of the interaction to be r e g = 98 ao and 
the eigenenergies using an energy-dependent pseudopo- 
tential [14| (dashed line in figure |3J). Both models agree 
to within a few percent, which is small compared to ex- 
perimental uncertainties. Further improvements taking 
into account more details of the atom-atom interaction 
in a two-channel model have been suggested 0, 0] and 
could be tested with our data. 

From a quantitative analysis of the spectra we obtain 
information about the occupancy of our lattice. We mea- 
sure the ratio between the atomic and the molecular 
peak heights in the spectra of the | — 5/2) atoms and 
determine the fraction of atoms that where bound in a 
molecule. Figure El shows the measured data for a lattice 
with a potential depth of 15 E r . The detected molecular 
fraction decreases for large values of 1/a, i.e. for deeply 
bound molecules, because of the small overlap of the ini- 
tial molecular and the final atomic wave function in the 
rf transition. From the analytical wave function of the 
molecular state in relative coordinates ip m (r) [ljj we cal- 
culate the overlap integral with the harmonic oscillator 
ground state ipho( r ) = {' K a\ )~ z / i exp(— r 2 /2af lo ) which 
determines the relative strength of the molecular tran- 
sition assuming that the center-of-mass motion remains 
unaffected by the rf transition. We fit the overlap in- 
tegral to our experimental data and obtain the fraction 
of molecules in the lattice to be (43 ±5)% (solid line in 
figure 131) - 

The observed fraction of molecules is primarily deter- 
mined by the filling of the lattice. To study the rela- 
tion between the temperature and filling we have numer- 
ically calculated the density of states for noninteracting 
fermions in an optical lattice including the Gaussian con- 
fining potential due to the transverse envelope of the lat- 
tice lasers. In the low-tunnelling limit we find that the 
density of states approaches p(E) oc E u with v = 1/2 
independent of the lattice depth [l^. The fraction of 
doubly occupied lattice sites in a 50:50 spin mixture de- 
pends on v and T/Tp. This makes the molecule frac- 
tion a quantity ideally suited for thermometry. Assum- 
ing that the occupation probability per spin state and 
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molecules formed. This allows for thermometry in the 
lattice, which has previously been unaccessible. The frac- 
tion of created molecules gives direct access to the num- 
ber of doubly occupied lattice sites in a two-component 
Fermi gas. Therefore it could be employed to character- 
ize a Mott insulating phase where the double occupancy 
should be strongly reduced. 

We would like to thank R. Diener, H. G. Katzgraber, 
and M. Troyer for insightful discussions, and OLAQUI, 
SNF and SEP Information Sciences for funding. 
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FIG. 3: The fraction of molecules detected by rf-spectroscopy 
at a potential depth of 15 E r . For weakly bound molecules 
(a^o/a < — 1) the dissociation works well since the overlap 
between the molecular and the atomic wave function is large. 
For deeply bound molecules aho/a > the detected molecular 
fraction is suppressed due to the vanishing overlap between 
the wave functions. The solid line shows the theoretical ex- 
pectation for a constant molecular fraction of 43%. 



lattice site i is < n.j < 1, the probability of finding 
two atoms with different spin state on a lattice site is 
nf. We identify the molecule fraction with the mean 
value of nf over the whole lattice. From a comparison 
with a numerical calculation |l8| we can conclude that 
the temperature of the atoms in the optical lattice is at 
most T/Tf — 0.45 ± 0.03. A similar result for our ex- 
perimental parameters was computed in |17| . This value 
gives an upper limit to the temperature since it assumes 
adiabatic formation of molecules at all doubly occupied 
sites and a perfect 50:50 mixture of the spin states. Dur- 
ing the ramp across the Feshbach resonance the density 
distribution might slightly change as compared to the ini- 
tial noninteracting case, which limits the accuracy of the 
temperature determination. 

We measure the lifetime of the molecules in the lattice 
to be on the same order of magnitude as of molecules in 
an optical dipole trap |3jJ. This is probably related to 
the tunnelling rate, which for our lattice parameters is 
on the order of 100 Hz and therefore the rate of three- 
body collisions is comparable to the optical trap. We 
expect that the lifetime of the molecules will significantly 
increase if the lattice depth becomes larger than 2>QE r . 
Another possible mechanism which could lead to a loss 
of molecules is photo-association induced by the lattice 
laser beams. 

In conclusion, we have studied molecules in a three- 
dimensional optical lattice. We have measured the bind- 
ing energies and find good agreement with the funda- 
mental theoretical model of two interacting particles in 
a harmonic potential well. Moreover, we have measured 
the filling of the lattice by determining the fraction of 
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